Ponderomotive forces and wave dispersion: two sides of the same coin 
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Presented here is a general view on adiabatic and resonant wave-particle interactions leading 
to a uniform description of nonlinear ponderomotive effects in very different environments, from 
low-temperature plasmas to relativistic plasmas or even atoms in laser light. Treating the wave- 
particle interaction as a classical mode-coupling problem, this theory shows the inherent connection 
between the ponderomotive forces and the properties of waves causing those forces. The adiabatic 
Lagrangians are derived for single particles and nonlinear waves, possibly carrying trapped particles, 
and yield both the dynamic equations and the nonlinear dispersion relations in the general case. 



1. Introduction. — A particle traveling in a rapidly 
oscillating wave field attains an average interaction en- 
ergy U acting as an effective "ponderomotive" potential. 
In low-temperature plasmas, U easily can be made com- 
parable or larger than the particle kinetic energy, thus 
producing barriers, which are generally asymmetric, ir- 
reversible, and exhibit quantumlike properties [1|. It was 
found recently that certain features of these barriers are 
akin to wave-particle interactions in also very different 
environments, from a channeling in fusion plasmas to 
laser traps for atomic cooling [l|-y and even relativis- 
tic laboratory and cosmological plasmas [5|-l7|. Elab- 
orating on these analogies could help understand the 
key underlying physics and also facilitate sharing tech- 
niques of wave and particle manipulation between differ- 
ent regimes. Thus, a unifying theory of ponderomotive 
interactions is needed. 

In Refs. [4. Isl. IsMlOj. a suitable analytical framework 
was developed by approaching wave-particle interactions 
as a variety of the classical mode-coupling problem. In 
this tutorial we highlight most general aspects of this 
framework, leaving out literature review and discussions 
on effects specific to particular physical systems. For 
those, one is referred to our earlier publications, and more 
explanatory papers will also follow. 

2. Single particle in a given field. — Let us start out 
with considering a particle in a given wave field, assum- 
ing that collisions and dissipation can be neglected. (If 
necessary, those can be included later as perturbations 
[ill .) Then the particle trajectory is found from the least- 
action principle: among trajectories x(i) starting at Xi 
at time ti and ending at X2 at time ^2 , realized is the one 
on which the action J^ Ldt is minimal. Here L is the 
particle Lagrangian, which we represent as L = {L) + L^, 
with (L) being the local time-average over the particle os- 
cillations, and L^ being the quiver part. On time scales 
large compared to the oscillation time scales, it is only 
(L) that contributes to the integral; hence, A = (L) acts 
as the slow-motion Lagrangian. 

To make use of this theorem, one needs to decompose 
the particle coordinate and velocity (x, v) into the slow 
variables (X, V) and the quiver variables (x, v), express 
the latter in terms of (X,V,t) [12'|, substitute the re- 



sult into L(x,v, i), and perform the time-averaging to 
obtain A(X,V,i), the time dependence retained being 
slow. Then, the motion equations are yielded in the form 
of Euler-Lagrange equations as usual; see, e.g., Ref. [13| . 
In general, the particle itself is an oscillator, i.e., has 
degrees of freedom additional to that of the translational 
"oscillation center" (OC) motion (x,v). (Examples here 
are gyromotion, intramolecular vibrations, etc.) Hence, 
it can be assigned one or more internal frequencies U = (p, 
with if being the corresponding canonical phase. If these 
oscillations are fast enough, the explicit dependence on 
If in the Lagrangian vanishes at time-averaging, so A can 
be written as a function of (x, v, fi, t); cf. Ref. [J]. Then, 
the canonical action J = i9aA is an adiabatic invariant, 
so one can reformulate the variational principle such that 
the internal degrees of freedom are omitted. Namely [5| , 
among trajectories x(i) starting at xi at time ti and 
ending at X2 at time t2, with arbitrary (pi and Lpi, realized 
is the one on which the action J^^ Cdt is minimal. Here, 
the reduced Lagrangian C, or the Routhian [ij] , reads as 



C^A-jn, f7 = f7(J,x,v,t), 



(1) 



so it is a function of (x, v, t) and depends on J paramet- 
rically. Then, the OC Hamiltonian equals 



H = P-Sr-C, 



(2) 



where P = c^v-C is the canonical momentum J5|. In par- 
ticular, dtH — dtH, so, in a stationary wave, T-L is an 
adiabatic invariant too [3|]. 

Notice that C, P, and Ti, can be of non-Newtonian form 
already in the nonrelativistic limit, thus giving rise to a 
plethora of unusual effects, e.g., the negative-mass effect 
[allOillSl- Still, the same formalism is also applicable in 
fully relativistic calculations [5|, |l3|, ll6[ . 

3. Invariants at resonant interactions. — Suppose 
now that the particle interacts resonantly with a station- 
ary wave, thereby contributing a periodic perturbation 
T-Lr^ to T-L. Then J and % are no longer conserved sepa- 
rately; yet, another general invariant can be found. To 
see this, let us extend the OC phase space by formally 
adding another degree of freedom {9, 1) corresponding 
to the oscillations at the wave frequency w. Hence the 



explicit time dependence in "Hr^ is replaced with the de- 
pendence on the phase 9, and for the action variable one 
obtains / = —T-L/u) [3[. Now our system formally con- 
tains two oscillators, (fi, J) and (w, /), so the generalized 
Manley-Rowe theorem applies [l7[. Specifically, suppose 
the resonance condition in the form, say, n^u ~ nQfl, 
where n^^ and riQ are some integers. Then, although 
J and / are no longer conserved separately, they are 
yet constrained to the so-called diffusion path [3|, llSl j. 
dl/n^ = -dJ/nn, yielding 

'H - {nu/nii) ujj — const. (3) 

Equation (jS]) generalizes the integral found in Refs. [19|, 
|20[ for asymmetric barriers in low-temperature plasmas 
and is also key to the cooling mechanism proposed in 
Refs. pj, l2|- The latter is similar to the celebrated Sisy- 
phus cooling [21|, showing inherent connections of the 
present theory with the physics of laser-atom interactions 
y. Yet, remarkably, the same conservation law [3| also 
facilitates a channeling in fusion plasmas 22| and thus 
suggests that akin techniques can be used to manipulate 
particles in very different environments. 

4. Ponderomotive forces from dispersion. — The 
angle-action representation for the wave variables turns 
out to be convenient also for calculating the pondero- 
motive Hamiltonian T-l itself, even when the particle dy- 
namics is considered in a prescribed classical field. For 
example, consider the wave as a mode with a prescribed 
spatial structure and suppose, for simplicity, that it is 
coupled adiabatically to a single particle. Then, J and / 
are conserved (the latter henceforth defined as the wave 
true action, unlike in Sec. 3), so what is affected by the 
interaction are the eigenfrequencies of the wave-particle 
system. In Ref. |8J , we showed that the shifts of these 
frequencies satisfy 



Sn = (9jt/)v, Slu = (diU)^, 



(4) 



where U = C — Cq is the interaction Lagrangian, and 
Co is the particle OC Lagrangian in the absence of the 
wave. (Particularly notice that the derivatives are taken 
at fixed V.) This means that knowing the wave micro- 
scopic frequency shift Suj per particle, one can, in princi- 
ple, find U. In other words, knowing the wave dispersion 
is sufficient to derive the ponderomotive force. 

In the simplest case when the interaction is linear, Slo is 
independent of /, so one obtains J7 w $, where ^ = 1 6u! 
is called the dipole ponderomotive potential. Since / per 
unit volume equals S/uj, where £ is the wave energy den- 
sity, one thereby obtains that, for species s, the pondero- 
motive potential $s is proportional to the derivative of 
w with respect to the particle density ris'. 



$., 



£ duj 
u dus 



(5) 



Equation ([S]) proves useful, e.g., for calculating the 
effect of ionization and recombination on homogeneous 



waves in cold plasma [23[. Yet, it can also be used to 
obtain $ in terms of the local field variables rather than 
those of an actual mode. Namely, substituting £ and the 
general dispersion relation for uj, one obtains [8] 



$ 



.E*-a-E, 



(6) 



where E is the electric field complex amplitude, and a. is 
the particle linear polarizability. (Alternatively, Eq. © 
is obtained by assigning to the particle a dipole moment 
q; • E, which characterizes the charge displacement from 
the OC trajectory [J|.) This also yields that 



•H«-Ho + $, 



(7) 



where "Ho is the OC Hamiltonian in the absence of the 
wave. In other words, the second-order (in E) interaction 
Hamiltonian is determined by the linear polarizability, 
the statement being known as the K-x theorem. In this 
approximation, the wave-particle energy is also separated 
automatically into the wave energy and the OC energy, 
with zero coupling between the two '8l|. 

If a particle is an oscillator, a exhibits a singularity 
near the resonance. Namely, at a; sa f7, one obtains 
<& OC |i?p/(a; — ri), always with a positive coefficient [4]. A 
nonlinear oscillator will yield a hysteretic U , as one finds 
similarly by integrating Eq. (U) for a model L [J]. Fi- 
nally, if multiple internal modes are present, coupling at 
beat resonances yields U matching the effective quantum 
potential seen by a two-level atom in a laser field. For 
example, in the simplest case of large enough detuning 
IS. — UJ — (fii — O2), the latter reads as 



U 






"1), 



(8) 



where fJ^ is the Rabi frequency, and ni are the corre- 
sponding occupation numbers [J]. 

Remarkably, these results are not limited to specific 
physical systems and ensure that ponderomotive forces 
in the absence of dissipation are of manifestly Lagrangian 
or Hamiltonian form. Hence, such forces conserve phase 
space automatically, unlike when they are derived ad hoc 
from time-averaging of the motion equations. 

5. Dispersion from, ponderomotive forces. — One can 
also invert the argument and derive uj from U instead. 
This can be done by summing contributions of individ- 
ual particles in the right-hand side of Eq. (|3]), so a differ- 
ential equation for uj is yielded. Alternatively, the wave 
dispersion is obtained in a nondifferential form as fol- 
lows. Along the lines of Sec. 2, the wave Lagrangian 
(Routhian) density is derived in the geometrical-optics 
approximation [9|, |24| , reading as 



\^cni) 






ns{Hs)f. 



(9) 



Here (£cm) = {E^ — B^)/{8t:) is the time-averaged La- 
grangian density of the electromagnetic field, with elec- 
tric and magnetic fields E and B including both wave 



and quasistatic fields, if any; summation is taken over 
distinct species s, and (Hs)/ are the corresponding OC 
Hamiltonians averaged over the local distributions fs ■ 

Understanding the wave properties hereby becomes 
straightforward. Namely, to obtain £, one needs to cal- 
culate only the single-particle Hamiltonians Us. The lat- 
ter is done along the lines of Sec. 2, thus rendering the 
Maxwell- Vlasov system unnecessary in many cases of in- 
terest, unlike in other existing theories [9|. In the sim- 
plest case, one can use the ponderomotive Hamiltonian 
in the dipole limit [Eqs. (jH]) and ([7])]; then. 



£ = £n 



1 
16^ 



E* -e-E- IBI' 



(10) 



with £o independent of the field variables. (Here we 
used that the linear dielectric tensor e equals 1 -I- 
47r^^ 715(0:5)^.) Hence, nondissipative linear waves are 
recovered immediately [Ml- Yet, Eq. © also describes 
nonlinear waves, nonperturbatively in the field amplitude 
a. The general dispersion relation is obtained by varying 
£ with respect to a(x, t); specifically. 



da£ = 0. 



(11) 



The envelope equation, or the action conservation the- 
orem 25| can be derived as well, by varying £ with re- 



spect to the wave phase 5(x,i). In the simplest case, 
Whitham's equation is recovered [26|: 



dt{d^£)-V-idk2) = 0. 



(12) 



However, we found that it can also be altered, namely, 
in the presence of particles trapped by the wave [2J] . A 
peculiar feature of such particles is that their OC loca- 
tions X are determined by the motion of the wave nodes 
and thus the trapped-particle density is connected with 
^. Also, the wave Lagrangian density £ attains a term 
that, albeit being a function of w and k, is independent of 
a. This special property of trapped-particle waves results 
in effects that r nay not be captured correctly by other ex- 
isting theories [2J|. More explanatory papers will follow 
that will describe our specific findings in this area. 

6. Summary. — Presented here is a general view on 
adiabatic and resonant wave-particle interactions lead- 
ing to a uniform description of nonlinear ponderomo- 
tive effects in very different environments, from low- 
temperature plasmas to relativistic plasmas or even 
atoms in laser light. Treating the wave-particle inter- 
action as a classical mode-coupling problem, this the- 
ory shows the inherent connection between the pondero- 
motive forces and the properties of waves causing these 
forces. Both are attributed to the same effect, namely, 
the eigenfrequency shifts in the wave-particle system. In 
particular, it is made clear how nonlinear ponderomo- 
tive energy in the dipole limit is responsible for the linear 
dispersion, and vice versa. Yet, the same formulation re- 
solves also essentially nonlinear effects, e.g., higher-order 



ponderomotive forces and effects of trapped particles on 
the wave properties in plasma. 
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